

full_universe =
  analysis_df %>%
  filter(final_score < 100) %>%
  mutate(fill_color = case_when(abs(forcing_f1) < 1 ~ "sample",
                                TRUE ~ "not_sample")) %>%
  ggplot(aes(x=final_score, y = forcing_f1, color = fill_color, alpha = fill_color, fill = fill_color)) +
  geom_point() +
  scale_color_manual(values = c("lightgrey", "black")) +
  scale_alpha_manual(values = c(0.4, 1)) +
  theme_bw() +
  theme(panel.grid.minor = element_blank(), 
        panel.grid.major.x = element_blank(),
        axis.line = element_line(colour = "black"),
        panel.border = element_blank(),
        legend.title = element_blank(),
        legend.position = "none") +
  ylim(-50, 50) +
  geom_hline(yintercept = 0, linetype = "dashed", color = "black") +
  ylab("Integrated Score (Difference)") +
  xlab("Integrated Score (Absolute)") +
  theme(text=element_text(size=12, family="lmroman"))
  

full_sample =
  analysis_df_subset %>%
  filter(final_score < 100) %>%
  mutate(fill_color = case_when(abs(forcing_f1) < 1 ~ "sample",
                                TRUE ~ "not_sample")) %>%
  ggplot(aes(x=final_score, y = forcing_f1, color = fill_color, alpha = fill_color, fill = fill_color)) +
  geom_point() +
  scale_color_manual(values = c("lightgrey", "black")) +
  scale_alpha_manual(values = c(0.4, 1)) +
  theme_bw() +
  theme(panel.grid.minor = element_blank(), 
        panel.grid.major.x = element_blank(),
        axis.line = element_line(colour = "black"),
        panel.border = element_blank(),
        legend.title = element_blank(),
        legend.position = "none") +
  ylim(-50, 50) +
  geom_hline(yintercept = 0, linetype = "dashed", color = "black") +
  ylab("Integrated Score (Difference)") +
  xlab("Integrated Score (Absolute)") +
  theme(text=element_text(size=12, family="lmroman"))

ggsave(filename = "./BKN survey/apsr_kuipers_replication_file/_4_outputs/figures/3_figure_a1_a.png", plot = full_universe, width = 6, height = 4)
ggsave(filename = "./BKN survey/apsr_kuipers_replication_file/_4_outputs/figures/3_figure_a1_b.png", plot = full_sample, width = 6, height = 4)

